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We study the Kondo Lattice Model (KLM) on a square lattice through a Hartree-Fock approxi¬ 
mation in which the local spins are treated semi-classically, in the sense that their average values are 
modulated by a magnetic wavevector Q while they couple with the conduction electrons through 
fermion operators. In this way, we obtain a ground state phase diagram in which spiral magnetic 
phases (in which the wavevector depends on the coupling constants and on the density) interpo¬ 
late between the low-density ferromagnetic phase and the antiferromagnetic phase at half filling; 
within small regions of the phase diagram commensurate magnetic phases can coexist with Kondo 
screening. We have also obtained ‘Doniach-like’ diagrams, showing the effect of temperature on the 
ground state phases, and established that for some ranges of the model parameters (the exchange 
conpling and conduction electron density) the magnetic wavevector changes with temperature, either 
continuously or abruptly (e.g., from spiral to ferromagnetic). 

PACS numbers: 71.27.-l-a, 71.10.Fd, 75.10.-b, 75.30.Mb 


I. INTRODUCTION 

The Kondo Lattice Model (KLM) [U [5] describes a 
system consisting of magnetic moments occupying each 
site of a regular lattice, and interacting with conduction 
electrons through a local exchange interaction, J; see, 
e.g., Refs.El andm The exchange coupling leads to two 
types of effects: the Kondo screening and the Ruderman- 
Kittel-Kasuya-Yosida interaction (RKKY). The former 
creates a hybridization between conduction electrons and 
local moments, which favours a paramagnetic Fermi liq¬ 
uid phase in which the local moments contribute to the 
Fermi surface; the latter creates an indirect interaction 
between the local moments, through the polarization of 
the conduction electrons, thus favouring a magnetically 
ordered state with a small Fermi surface accommodating 
solely conduction electrons. As first discussed by Do- 
niach [T], these two effects can coexist, hence compete 
with each other, leading to a quantum phase transition 
from a magnetically ordered state to a paramagnetic one. 

This competition between Kondo screening and mag¬ 
netism, including the existence of a quantum critical 
point, is known to occur in several heavy-fermion ma¬ 
terials [1]; in view of this, it is generally believed that 
the KLM provides an adequate description of some as¬ 
pects of these materials. Further, early mean-field ap¬ 
proaches to the Kondo-lattice Hamiltonian mm already 
pointed out that for small screening, J < W {W is the 
bandwidth), the magnetic ground state is either ferro¬ 
magnetic or antiferromagnetic, respectively for electronic 
densities 0 < Uc < 0.6 and 0.6 ric < 1. This im¬ 
mediately connects with an interesting class of materi¬ 
als, the borocarbide family [5], which has the chemical 
composition i?r 2 B 2 C, where R represents a rare-earth 
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element and T is a transition metal. The rare earth 
usually contributes with one local magnetic moment per 
unit cell, which, in turn, interacts with the conduction 
electrons. Some members of the family display coexis¬ 
tence between superconductivity and magnetism EHin], 
the interplay of which is a subject of current interest as a 
unifying link with the iron pnictides m, and, possibly, 
with the cuprates [laila]. Setting aside the supercon¬ 
ducting behaviour of some borocarbides, one notes that a 
wide variety of magnetic orderings (or modes) have been 
found: depending on the particular combination of R 
and T, one finds ferromagnetism, antiferromagnetism, as 
well as spin-density waves, and multiple Q-wavevectors, 
commensurate and incommensurate IS m]. This vari¬ 
ety of magnetic modes has been further scrutinized with 
the synthesis of borocarbides with variable proportions 
of transition metals, as in Tb(Coa;Nii_a;) 2 B 2 C [Refs. US 
andfTH] and Ho(Coa;Nii_a;) 2 B 2 C [H ef.fTT] . Notwithstand¬ 
ing the fact that the KLM does not incorporate explic¬ 
itly some aspects of the borocarbides (such as crystal 
field effects), one may wonder whether it can be used as 
an effective model to describe the evolution of magnetic 
modes with the band filling. 

At any rate, the KLM is an interesting model in its 
own right, and a great deal of theoretical effort has been 
invested to unveil its overall properties. One-dimensional 
topologies are amenable to unbiased methods such as the 
Density Matrix Renormalization Group (DMRG), and, 
indeed, ferromagnetism and spiral magnetic phases can 
be stabilized in the linear chain, in different regions of 
the ground state parameter space (J, Uc), where ric is 
the density of conduction electrons HU; further, the two- 
leg version of the model exhibits quasi-long-range mag¬ 
netic order, with the magnetic wavevector displaying a 
well defined dependence with the electronic density m- 
In two dimensions, auxiliary-field Quantum Monte Carlo 
(QMC) results are available at half filling 123 El]; away 
from half filling, variational QMC has been used, but only 
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antiferromagnetic solutions were probed |22l I23j . In view 
of the intrinsic difficulties of those methods to extract 
comprehensive and simultaneous information (i) about 
the various magnetic arrangements, (ii) about the inter¬ 
play between Kondo screening and magnetism, (iii) in 
two- and three dimensions, (iv) for all conduction elec¬ 
tronic densities (from 0 to half filling), and (v) the ef¬ 
fects of temperature, many mean-field implementations 
have been used over the years, tackling some of these is¬ 
sues. While several studies only allowed for para-, ferro- 
and antiferromagnetic phases [113 HMH], the possibil¬ 
ity of Q-dependent magnetic modes was considered in 
Ref.[23 in which case the local spins were treated clas¬ 
sically, thus precluding the analysis of the coexistence 
of Kondo screening with magnetic phases. By contrast, 
this coexistence has so far been examined only in con¬ 
junction with ferromagnetic (FM) or antiferromagnetic 
(AFM) phases [23H26] . 

Therefore, a mean-field investigation of the KLM tak¬ 
ing into account both generic magnetic orderings and 
the effect of Kondo screening is clearly in order. With 
this in mind, here we use a ‘semi-classical’ approach, in 
which the local moments display a Q-dependent average 
magnetization, while they are also expressed in terms of 
fermionic operators allowing us to define a hybridization 
‘order parameter’ as a measure of the Kondo screening. 
In this quest, we are led to minimize the free energy also 
with respect to the magnetic wavevector, in order to es¬ 
tablish the dependence of the stable Q values with J, 
Uc-, and the temperature, T. Since this is more readily 
carried out in two dimensions than in three, we choose to 
consider here the KLM on a square lattice. Our main re¬ 
sults can be summarized in the form of a phase diagram 
with many magnetically-ordered phases which evolve into 
a screening-dominated (Kondo) one, but going through 
intermediate regions in which they coexist. 

This paper is organized as follows: The model is pre¬ 
sented in Sec.|n| together with highlights of the standard 
Hartree-Fock approximation, the details of which can be 
found in the Appendix. In Section III we discuss the 
results for the ground state, while Section jlV] is devoted 


to finite temperature behaviour. And, finally, Section |V| 
summarizes our findings. 


II. MODEL AND METHOD 


The Kondo lattice model is described by the Hamilto¬ 
nian 


^ E + H-C-) + J 5] S, • s^, (1) 


where the sums run over sites of a two-dimensional square 
lattice, with (z, j) denoting nearest-neighbor sites. The 
first term represents the hopping of conduction electrons, 
where (czo-) is the creation (annihilation) operator for 
an electron on site i with spin cr, and H.c. stands for 
hermitian conjugate of the previous expression; t sets the 
energy scale. The second term represents an interaction 
between local moments and conduction electrons, where 
J > 0 is the coupling strength, and and s)’ are the spin 
operators for the local moment and conduction electrons, 
respectively. 

In order to set up a Hartree-Fock approximation, we 
write the spin operators in a fermionic basis as 

= J E (2) 

a./3=± 

and 

= b E (3) 

Q!,/3=± 

with cra,i 3 denoting Pauli matrix elements, and (fia) 

being the creation (annihilation) operator for a localized 
electron with spin a on site i. Following the procedure 
outlined in the Appendix, the Hartree-Fock Hamiltonian 
becomes 


nMF = -tJ2 (els- + H-C.) + ^E+ tE(V“-(V/') + 


3J 

Y 


- ¥ E - A E(V") ■ (E) + ^ E(E) (E/) - ^ E(s*) ■’ (4) 


with the definitions 


v« = v<>; 


a,/3=± 


where I is the identity matrix, and 

= = i (6) 

a,/3=± 

Following the nomenclature introduced in Ref. l29L we re- 
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fer to and as singlet hybridization operator, and 
to Vic and 'Vif as triplet hybridization operators. 

In order to analyse the stability of planar spiral mag¬ 
netic phases, the mean value (S^) is taken as classical. 


where /? = fcs = 1 throughout this paper. The 

effective fields Vb, m°, /i, e/, and Q are to be de¬ 

termined self-consistently by minimizing the Helmholtz 
free energy 


(Si) = [cos(Q-Ri) ,sin(Q-Rj) ,0], (7) 


with 


Q — (fe) Qy) (8) 

being the magnetic wavevector, and Ri the vector po¬ 
sition of site i on the lattice. By the same token, we 
choose 


(Si) = -mc[cos(Q-R*) ,sin(Q-Ri) ,0], (9) 

where the minus sign above reflects the local antiferro¬ 
magnetic coupling between the local moments and the 
conduction electrons. 

The singlet hybridization terms can be taken as 

(^.c) = (Kj) = -Vo, (10) 

and the mean values of the triplet hybridization operators 
are similarly assumed to be given by 

{V^c) = (V,/) = Vo'[cos(Q-R,),sin(Q.R,),0]. (11) 


The electronic density, and the number of local 
moments per site, respectively expressed by 


^ t _ 

i(7 


( 12 ) 


and 


- V f =1 

^ j laJ 1(7 ’ 


(13) 


are imposed as constraints through the method of La¬ 
grange multipliers. The latter constraint is enforced on 
average, which seems to be unavoidable in mean-field 
treatments; as pointed out in Ref.|22l this may restrict 
analyses on the character of the Fermi surface, as far as 
being hole-like or electron-like, large or small. However, 
our main purpose here is to gain insight into the stabi¬ 
lization of different magnetic modes, so that the tradeoff 
justifies imposing the constraint in its weaker form. 

As discussed in the Appendix, after substituting 
Eqs. ([^-([TT|) in the mean-field Hamiltonian, Eq. Q, im¬ 
posing periodic boundary conditions, and performing a 
discrete Fourier transform, we obtain our working Hamil¬ 
tonian, Eq. (A.17). It is represented by a 4 x 4 matrix, 


which can be straightforwardly diagonalized, leading to 
the bands {n = 1,..., 4). 

The Helmholtz free energy then becomes 


F = ln(l-be-'^^^ 

n,k 


const. 


(14) 



The resulting nonlinear coupled equations are solved nu¬ 
merically, using standard library routine packages, with 
the aid of the Hellmann-Feynman theorem. 

At this point some comments are in order. First, 
we should mention that we have tried to include addi¬ 
tional modulations to the field amplitudes, e.g., |(Si)| —>• 
m/(Ri) = -I- TOj cos (Q-Ri), and similarly for |(s°)|, 

|(V)| and |(V)|, but, as it turned out, the most sta¬ 
ble solution always yields = Vi = V( = 0 

for all ranges oi J/W and Uc considered. Second, at¬ 
tempts to consider different Q’s for any of (S^), (sj) and 
V amount to a much harder minimization procedure, and 
have led either to unphysical results, such as spin ampli¬ 
tudes larger than 1/2, or to trivial mean-field solutions. 


III. GROUND STATE BEHAVIOUR 


The set of nonlinear coupled equations, Eq. (15), is 
solved numerically for each pair of (ric, J), by fixing the 
temperature (T = 0, for the time being) and the elec¬ 
tronic density, and by letting the exchange coupling to 
vary. Figures to show the behaviour of the order 
parameter amplitudes, Eqs. 0-(II3 , as functions of the 
exchange coupling J (in units of the bandwidth, W = 8t) 
for different doping levels (i.e., nc < 1). The figures also 
display the behaviour of the magnetic wavevector with 
J/W: as we will see, the magnetic modes and the co¬ 
existence with Kondo screening depend strongly on the 
electronic density. 

At half filling, the system is known to be an insulator 
for all J/W, but a quantum phase transition between an 
antiferromagnetic state and a spin singlet takes place at 
{J/W)c — 0.4. Care must be taken when comparing this 
estimate with those of Refs. uniisii and [Ml since their 
working Hamiltonians (i.e., after some decouplings or ef¬ 
fective Hamiltonians are introduced) is somewhat differ¬ 
ent from ours, Eq. (A.17). Nonetheless, a rough corre¬ 
spondence can be worked out from which the exchange 
coupling J used in those works is one half of the one we 
use here; with this proviso, our estimate for {J/W)c is in 
good agreement with those of Refs.jMHHl andIMl 

Figure [^a) shows the results for Uc = 0.90. In the 
weak-coupling regime, there is no hybridization (hence 
no Kondo screening); the local-moment amplitude is not 
affected by the exchange, while the amplitude of the con¬ 
duction electron magnetization increases steadily with 
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FIG. 1. (Colour online) (a) Order parameter amplitndes as 
fnnctions of the Kondo exchange conpling J (in units of the 
bandwidth, W = 8t): local moment, nif (fnll black line), 
conduction electron magnetization, (dashed blue line), 
singlet hybridization, Vo (dashed-dotted red line) and triplet 
hybridization, Vq (full orange line), (b) Magnetic wavevector 
component g as a function of JjW', the spiral magnetic phase 
is described by the wavevector Q = (rr,?). All data are for 
conduction electron density = 0.90, and zero temperature. 


J/W. The stable magnetically ordered phase corre¬ 
sponds to a spiral arrangement with wavevector Q = 
(tTjq) [or, by symmetry, with (q,7r)], and from Fig. 0b). 
we see that q first decreases slightly with J/W, hence 
drifting away from the Neel case; that is, the increase 
in the Kondo coupling by itself cannot drive the sys¬ 
tem into an antiferromagnetic state. However, when 
J/W « 0.35, hybridization abruptly sets in, causing a 
sudden decrease in both magnetic amplitudes; the mag¬ 
netic mode also changes abruptly, stabilizing an AFM 
phase, with Q = (tt, tt), which coexists and competes 
with Kondo screening. Further increase in J/W en¬ 
hances the singlet hybridization which, in turn, steadily 
suppresses the magnetic amplitudes, both vanishing at 
J/W « 0.375. The triplet hybridization amplitude Kq 
also vanishes at J/W « 0.375, tracking the suppression 
of the magnetization. Beyond this point, there is only 
a paramagnetic phase with non-zero Vq, usually referred 
to as the Kondo phase. One can also see from Fig.[^a) 
that the order parameters are discontinuous across the 
lower transition (into the coexistence region), and con¬ 
tinuous at the second transition, into the screened-only 
(Kondo) phase. The most stable ground state therefore 
corresponds to V/ <^Vo, which indicates that the compe¬ 
tition with magnetism is almost entirely due to the singlet 
hybridization; accordingly, from now on hybridization ef¬ 
fects will only be associated with Vq- 



FIG. 2. (Colour online) Same as Fig.j^ but for ric = 0.60. 

Moving on to tic = 0.60, we see from Fig.j^a) that 
the order parameters behave in a way similar to the case 
with ric = 0.90, including the order of the transitions; 
in addition, a similar magnetic mode with Q = (tt, q) is 
stabilized in this case. As shown in Fig.j^b), q also de¬ 
creases with J/W in the unscreened region, though with 
the important difference that in the coexistence region it 
is the (tt, 0) mode which dominates. Later on, we will 
discuss the behaviour of Q as a function of ric, for hxed 
J/W. 



J/W 


FIG. 3. (Colour online) Same as Fig.[^ but for = 0.35. 
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FIG. 4. (Color online) Same as Fig.[^ but for = 0.25, and 
now the magnetic wavevector is Q = (0,g). 
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FIG. 5. (Color online) The ground state phase diagram, 
Kondo exchange versus electronic density. The magnetic 
phases are designated by their magnetic wavevectors, Q = 
{qx,<ly), and ‘Kondo’ denotes a phase in which screening, as 
measured by the hybridization, is present (see text). Solid 
and dashed lines respectively represent continuous and dis¬ 
continuous phase transitions. 


Figure shows the corresponding analysis for Uc — 
0.35. We see that the stable magnetic phase now has 
a wavevector Q = (tt, 0) in the unscreened region. In 
the coexistence region, the wavevector first stabilizes in 
a mode Q = (<7,0), before becoming FM, Q = (0,0), 
as showed in Figure [^(b); see also Fig.[^ To the best 
of our knowledge, this is the hrst time that the coexis¬ 
tence of a spiral incommensurate magnetic mode with 
the Kondo phase is predicted within a static mean-field 
analysis; more on this coexistence later. 

Further decrease in the electronic density, e.g., for ric = 
0.25, leads to a spiral magnetic phase with Q = {q, 0); see 
Fig.lU While in the unscreened phase one finds a mono- 
tonically decreasing q{J/W), in the coexistence region a 
uniform FM phase [i.e., one with Q = (0,0)] is stabilized. 
Interestingly, while the transition into coexistence (which 
occurs at J/W « 0.51) is still of first order, here we see 
that, unlike what we have discussed so far, the transition 
to pure Kondo behaviour, occurring at J/W « 1.02 is 
also discontinuous. 

Similar analyses were performed for other values of 
ric, the results of which are summarized in the phase 
diagram of Fig.[^ In line with previous mean-field ap¬ 
proaches lan], we see that at low densities a saturated 
ferromagnetic phase is stable, while at half filling it is 
an antiferromagnetic phase which is the stable one. On 
the other hand, we have established that the evolution 
of magnetic modes with the electronic density is much 
smoother than hitherto assumed; the diagram of Fig.[^ 
also shows that magnetically ordered phases can still 
withstand some screening, though no trace of magnetism 
is found deep in the Kondo phase, as expected. 

We now discuss these aspects in turn, starting with the 


magnetic ordering in the absence of Kondo screening. As 
the electronic density increases from zero, a spiral mod¬ 
ulation develops in one of the lattice directions, say, the 
cc-direction, while the same modulation is repeated along 
the y-direction, for which qy = 0. Figure [^illustrates the 
evolution of the modulation vector with^e density, for 
fixed J/W: the modulation along x becomes staggered 
{<lx = and so remains, even as the density increases 
slightly up to Tic 0.4. As ric continues to increase, the 
modulation along y starts changing until it also reaches 
qy = n close to half hlling. 



FIG. 6. (Color online) The magnetic wavevector Q = {qx,qy) 
as a function of electron density, for fixed J/W = 0.125. 
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FIG. 7. (Color online) Internal energy contour map as a 
function of the magnetic wavevector, for ric = 0.60 and 
J/W = 0.175. 



q In 


FIG. 8. (Color online) Polarization contour map as a function 
of magnetic wavevector, for the same parameters as in Fig.[^ 


The presence of spiral magnetic phases has been found 
for the one-dimensional KLM [18], in addition to ferro¬ 
magnetic ‘island’ phases; while our results reveal that 
spiral phases also occur for the square lattice, no ferro¬ 
magnetic island states led to minimal free energies in the 
present case. The continuous change in one of the com¬ 
ponents of the magnetic wavevector can be attributed to 
a distribution of conduction electrons preferentially along 
one of the lattice directions; similar effects have been ob¬ 
served in the DMRG study of the 2-leg Kondo ladder 
m-, that is, the Kondo lattice seems to develop a stripe 
structure. The correspondence with the 2-leg Kondo lad¬ 
der goes even farther: fits to the linear portions of q{nc) 
in Fig.yield q and q « 271^, respectively near half¬ 
filling, and near ric ~ 0.35. This should be compared 
with Q = (ric, l)7r and Q = (2nc,0)7r, for ric ^ 0.5 and 
n-c ^ 0.5 respectively, for the ladder m- On the other 
hand, there is a noticeable difference for a range of den¬ 
sities near ric ~ 0.4, where, for the square lattice, the 
commensurate phase (tt, 0) stabilizes; see Fig.[^ 

Let us now compare these predictions with 
the experimental data for the borocarbide alloy 
Tb(Coa;Nii_a;) 2 B 2 C [Refs.lTSl and [T^. The planar 
magnetic arrangement in TbNi 2 B 2 C is a spin-density 
wave, with a wavevector close to (7r,0) [RefjT^j, so that 
we can represent this compound by the ric ~ 0.4 point 
in Fig.[^ Assuming the primary effect of the gradual 
substitution of Ni by Co is a decrease in the number 
of conduction electrons (hence of ric), Fig.|^ correctly 
predicts that the alloy evolves towards a saturated 
ferromagnet in the opposite limit of 100% Co; the 
comparison cannot be made for intermediate dilutions, 
since partial replacement mainly affects the modulation 
along the c-axis. In the corresponding case of the Ho 
alloys, the planar arrangement is ferromagnetic for all Co 
concentrations, x, while the magnetic modulation along 
the c-axis is strongly dependent on x. m However, since 
Ho(CoxNii_ 2 ,) 2 B 2 C is superconducting below x = 0.03, 


one expects electron-electron interactions to play an 
important role in the ensuing magnetic arrangement, 
even when the system is not superconducting. While 
these effects are certainly absent in the simple model 
considered here, the capture of the trend observed in the 
Tb alloys may be taken as an indication that the KLM 
is a viable starting point to describe the magnetism in 
this class of materials. 

Further insight can be gained by discussing the rela¬ 
tive stability between AFM and FM phases. Figure 
shows a contour map of the internal energy A as a func¬ 
tion of magnetic wavevector {qx,qy), for J/W = 0.175 
and ric = 0.60. The map is obtained by minimizing the 
energy, Eq. (15), with respect to all variables, but qx and 
qy. From Fig.j^ we see that when Q = (7r,7r) and (0,0), 
the internal energy reaches its largest values, showing 
that for this choice of {ric, J/W) the most stable mag¬ 
netic arrangement in the ground state is neither AFM 
nor FM. The minimum of the internal energy actually oc¬ 
curs for {q, tt) [and, by symmetry for (tt, q), as well], with 
q/iT ^ 0.55, thus providing us with an explicit example 
showing that many different magnetic arrangements may 
be closely separated in energy. Figure]^ shows the corre¬ 
sponding contour map for the polarization, m/. We first 
note that the polarization is in opposite phase in relation 
to the internal energy. For this choice of {ric, J/W), the 
hybridization is zero (see Fig.j^, and the system is domi¬ 
nated by the RKKY interaction, so that the polarization 
is the sole driving force to magnetism. 

As far as the order of the density-driven transitions 
is concerned, we should add a few comments. First, we 
note that the transition is necessarily discontinuous if Q 
changes abruptly at the boundary; as shown in Fig. i Q 
does not suffer any discontinuity in the whole range of 
ric- Further evidence comes from Fig.j^ which shows the 
internal energy as a function of q, as the transition from 
{•7T,q) to (7r,0) is approached, for fixed J/W and varying 
ric'. we see that the two global minima go continuously 
to zero as the critical point is approached. Analogous 
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FIG. 9. Internal energy as a function of magnetic wavevec- 
tor as the zero-temperature transition from (tt, q) to (tt, 0 ) 
is approached by varying the electronic density, Uc, for fixed 
J/W = 0.125. The curves were shifted for easier comparison. 
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behaviour occurs for all other transitions at fixed JjW. 

We now discuss the coexistence of Kondo screening and 
magnetically ordered phases. First, we note that this co¬ 
existence is restricted to moderate degrees of screening, 
as indicated in Figs.to and summarized in Fig.[^ In 
line with other mean-field analyses [ 2 ni - l 2 H] . our results 
show that close to half filling, i.e., for ric ^ 0.86, Kondo 
screening coexists with an AFM mode, while for low elec¬ 
tronic densities (0.15 ric 0.31) this coexistence oc¬ 
curs with the FM configuration. However, our approach 
allows us to go further, and establish that for interme¬ 
diate electronic densities, 0.31 < < 0.68, coexistence 

is possible with phases other than FM and AFM; see 
Figs. ii and[^ A commensurate magnetic phase, with 
wavevector Q = (tTjO), stabilizes into the coexistence re¬ 
gion for the range 0.36 0.68. More interestingly, 

for the range 0.31 ’^ric 0.36 such coexistence is with a 
spiral incommensurate magnetic phase. 


Let us now focus on the coexistence between Kondo 
and spiral phases. Figure 10 (a) compares the internal 
energy as a function of q/Tr (for fixed JjW') when the 
triplet hybridisation amplitude (Kq) is constrained to be 
zero, and (b) when it is allowed to be non-zero. In the 
former case, the energy is minimum at the ferromagnetic 
mode Q = (0,0), while when this constraint is relaxed a 
mode with g/ 7 r « 0.25 becomes the most stable one. It 
should also be stressed that, by contrast, the final min¬ 
imisation outcome for the remaining auxiliary fields, m°. 


m° and Vq, is hardly affected by whether Vq is zero or 
non-zero. We conclude that the appearance of modes 
with (7 yf 0, TT in the region of coexistence between Kondo 
screening and magnetic order is directly related to a mod¬ 
ulation of the hybridisation with the magnetic wavevec¬ 
tor Q. Evidently, as J/W varies, the value of q which 
minimises the energy also varies; see Fig.[^(b). Another 
subtle aspect is that the resulting amplitude of the mod¬ 
ulated hybridisation is weak (typically V//Vq ~ 0.1), so 
that coexistence involving Kondo screening and modes 
with either g 0 or yf tt only occur in small portions of 


FIG. 10. Internal energy as a function of magnetic wavevector 
Q = {q, 0) at ric = 0.35 and for fixed J/W — 0.59, in the cases 
where (a) the triplet hybridisation term (Vq) is forced to be 
zero, and (b) when it is allowed to be non-zero, in the coexis¬ 
tence region. The thin black dashed line is the ferromagnetic 
internal energy for both cases. 


the diagram, the precise location of which would demand 
a much more elaborate analysis; suffices to say, for our 
purposes here, that coexistence with 9 yf 0 , 7 r is indeed 
possible. 

Finally, we examine the transition to the Kondo phase, 
which marks the disappearance of magnetism. As shown 
in Figs. BUI and the transitions to the pure Kondo 
regime from both the AFM-|-Kondo phase and from the 
(0, 7 r)-|-Kondo phase, are continuous. By contrast, the 
transition from FM-|-Kondo is discontinuous, as it can be 
seen from Fig. El in which we fix the electronic density 
as Tie = 0.30, and plot the internal energy as a function 



FIG. 11. Internal energy as a function of magnetization at 
ric = 0.30. The internal energy behaviour before {J/W = 
0.85), on {J/W ~ 0.866) and after {J/W = 0.881) the tran¬ 
sition point to pure Kondo are represented by black (solid), 
red (dashed) and blue (dash-dotted) lines, respectively. The 
curves were shifted for easier comparison. 
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FIG. 13. (Color online) Temperature dependence of the vari¬ 
able wavevector component q, for three distinct combinations 
of (ric, JjW). The magnetic wavevector in each panel is given 
by (a) Q = (q,0), (b) Q = (tt,?), and (c) Q = (vr,q). 


FIG. 12. (Color online) ‘Doniach diagram’ for = 0.25. 
Solid lines represent continuous transitions, while broken 
(dashed and dotted) lines denote discontinuous transitions; 
red lines mark the onset of hybridization, while black lines 
mark the disappearance of magnetic order. PM(N) stands 
for normal (i.e., no hybridization) paramagnetic phase. The 
dot in the boundary between FM-|-Kondo and Kondo phases 
marks the tricritical point; the insets show the behaviour of 
the Helmholtz free energy as a function of the average local 
moment, along the first-order boundary. 

of the local moment amplitude m°, for values of J/W 
near the transition point, {J/W)c « 0.866. Within the 
coexistence region, J/W = 0.85 < (J/W)c, the inter¬ 
nal energy displays global minima at « ±0.4, sig¬ 
nalling a FM state, together with a local minimum at 
= 0. At the transition point, this local minimum be¬ 
comes degenerate with those for which ~ ±0.4, and 
for J/W = 0.881 > {J/W)c, the minimum at vnPj = 0 
becomes the most stable one: the transition is therefore 
discontinuous. Similar discontinuous behaviour of the 
internal energy is found for the direct transitions (i.e., 
without going through coexistence regions) to the Kondo 
phase from the Q = (0,0) (for ric < 0.15) and Q = ( tt , q) 
(0.5 < nc < 0.85) phases; see Fig.[^ One should note 
that in Ref.|26|the transition FM±Kondo to Kondo was 
found to be continuous; possible sources for this discrep¬ 
ancy may lie in either the constant density of states used 
in that work, or to the fact that the mean-field implemen¬ 
tation differs from ours, especially for the FM solution 
(see the Appendix). 


IV. FINITE TEMPERATURES 

In analysing the behaviour at finite temperatures we 
fix ric, while the temperature and the exchange coupling 
are allowed to vary. Similarly to what we did for T = 
0 , we examine the temperature dependence of the order 
parameters to determine the phase boundaries; we also 


determine the temperature dependence of the magnetic 
wavevector. 

Figure [^summarizes our findings for Uc = 0.25 in the 
form of a ‘Doniach diagram’, in which several details on 
the nature of the magnetically ordered phases can now be 
unveiled. First, we note that for this electronic density, 
the unscreened magnetic phase is actually a spiral phase. 
More interestingly, the magnetic mode Q = {q, 0) is such 
that q displays a temperature dependence, as shown in 
Fig.j^a), for a fixed J/W = 0.56. In this case, the mag¬ 
netic mode is hardly dependent on the temperature in 
the unscreened region, but an abrupt change occurs as 
soon as the temperature drives the system into the coex¬ 
istence region, where, in this case, ferromagnetism sets in 
over a temperature interval; for a slightly smaller value, 
say J/W = 0.5, there is a noticeable temperature depen¬ 
dence of q with T (not shown), but the range of tem¬ 
peratures in which the FM phase exists is quite smaller 
than the one shown. The order of the transitions along 
the border between FM and PM Kondo phases changes 
from continuous (at higher temperatures) to discontinu¬ 
ous (lower temperatures); the insets show the evolution 
of the free energy along the first-order boundary, and at 
the tricritical point. For completeness, one should men¬ 
tion that a similar phase diagram was obtained in Ref. l3II 
for a single electronic density, ric = 0.2; however, since 
no spiral phases were considered there, the evolution of 
the wavevector Q with the temperature could not be es¬ 
tablished. 

The features brought about by the spiral phases also 
manifest themselves at other densities. Figure [l4| shows 
data for ric = 0.6, in which case the spiral phases in¬ 
volved at T = 0 are those for Q = {ir^q), with q de¬ 
creasing as J/W increases (see Fig.[^. One notes that 
the magnetic boundary between (7r,0)±Kondo and the 
Kondo phase is completely detached from the boundary 
between (7r,g) and PM(N); therefore, one can go from 
(tt, q) to the Kondo phase without an intervening coex¬ 
istence region, simply by raising the temperature. As 
Figure [I^b) shows, for J/W = 0.3 the effect of tern- 
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FIG. 14. (Color online) Same as Fig. 12 but for ric = 0.6. 


perature is to increase g, moving towards antiferromag¬ 
netism. On the other hand, the coexistence region which 
appears for 0.4 < J jW < 0.54 involves the ‘striped’ 
phase Q = (tt, 0). 

For ric = 0.9, the phase diagram is shown in Fig.flS] 
Unlike the previous cases, by increasing the temperature 
one can now smoothly interpolate from (Tr,^) to (7r,7r), 
provided the value of J jW lies entirely to the left of the 
first order line (red dashed curve in Fig. |15| ). By con¬ 
trast, if one chooses, say J jW = 0.28 as in Fig.[T^c), 
one obtains a reentrant behaviour for the AFM mode, 

Q = (7’‘,7r). 



J/W 


FIG. 15. (Color online) Same as Fig. |12| but for = 0.9. 


V. CONCLUSIONS 


In conclusion, we have analysed the Kondo Lattice 
model on a square lattice, using a semi-classical-spin ap¬ 
proach within a Hartree-Fock approximation. This al¬ 
lowed us to probe the presence of spiral magnetic modes, 
which, for some ranges of parameters, turned out to be 
more stable than the ferromagnetic, antiferromagnetic, 
and paramagnetic modes, usually assumed to be the only 
ones at play. The presence of spiral phases is in line with 
DMRG calculations for the one-dimensional case[THl, as 
well as for the two-leg ladder [19], so that one may expect 
they are not an artifact of the approximations employed 
here. Accordingly, we have obtained a ground state phase 
diagram in terms of the Kondo coupling, JjW {W is the 
bandwidth), and the conduction electron density, tlc- As 
ric varies from 0 to 1 (half-filling), the weak- to mod¬ 
erate coupling region displays a variety of incommensu¬ 
rate phases [i.e., with continuously changing magnetic 
wavevectors, in the form Q = (tt,?), ( 9 , 0 ), or ( 9 , 9 )], 
in which Kondo screening is absent. Recent dynami¬ 
cal mean-field theory (DMFT) calculations on the KLM 
away from half filling (see Ref.[32|) show incommensurate 
spin-density waves on both the small- and large Fermi 
surface regions; our results are in good agreement in the 
unscreened region. 

We have also found that magnetic phases can coex¬ 
ist with some degree of screening, and that the change 
in magnetic wavevector with the Kondo coupling occurs 
even in the region of coexistence with the Kondo phase. 
Such a coexistence with incommensurate magnetic modes 
occurs in just a tiny region of the phase diagram, be¬ 
ing related to the modulation of the triplet hybridisation 
term. It seems that a modulated hybridisation is an im¬ 
portant ingredient for the stabilisation of magnetic modes 
other than Q = (0,0) or Q = ( 7 r, 0 ) in the coexistence 
region. On the other hand, for sufficiently strong cou¬ 
pling, screening dominates and magnetism is suppressed. 
While completing this work we became aware of Ref. [33| , 
in which the decoupling of the Kondo term in singlet and 
triplet hybridisations is the same as ours, but, unlike our 
present framework, the magnetic wave vector Q was not 
left as a free parameter to be determined by minimisa¬ 
tion of the free energy. Since Li et al. [33] set it as (tt, tt), 
they could not discuss spiral phases; instead, they con¬ 
sidered the effects of an additional hopping term between 
next-nearest neighbours. 

We have also discussed the behaviour at finite tem¬ 
peratures in the form of ‘Doniach-like’ phase diagrams 
T X JjW, for hxed electronic densities. We have es¬ 
tablished that unscreened spiral magnetic phases can be 
found in the low-temperature and small-coupling por¬ 
tion of the phase diagram. Within each unscreened 
phase, the magnetic wavevector in general increases with 
temperature, until it reaches the hrst-order transition 
to the phase of screened magnetic order, when abrupt 
changes in 9 may occur; interestingly, the coexisting mag¬ 
netic mode is always commensurate, leading to antifer- 
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romagnetic, ferromagnetic, or striped phases. As ex¬ 
pected, at sufficiently high temperatures only param¬ 
agnetic phases survive, though they can be either un¬ 
screened or screened, depending on the magnitude of the 
Kondo coupling; however, as pointed out before [5T|. the 
sharp transition between these two regimes is expected to 
become a crossover if fluctuations were taken into account 
beyond a mean-field treatment. The results presented 
here suggest that through a judicious choice of parame¬ 
ters, Kondo lattice systems may allow for temperature- 
driven switching between detectable magnetic modes. 
Notwithstanding the fact that the present results have 
been obtained for a square lattice, we have found that the 
evolution of magnetic modes with the model parameters 
shares common trends with borocarbides family of mate¬ 
rials, so that the Kondo lattice model should provide an 
adequate description of their non-superconducting prop¬ 
erties. 

We close with a cautionary remark. The prediction of 
ordered states at finite temperatures in two-dimensional 
systems with continuous symmetry is certainly a draw¬ 
back of mean-field approximations. Nonetheless, ground 
state features like the continuous variation of Q with 
the electronic density (as described in the simple pic¬ 
ture above) should be present in both two- and three¬ 


dimensional systems; our results also broadly suggest 
how the temperature would influence the magnetic modes 
in three dimensions (or weakly-coupled two-dimensional 
layers). 
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Appendix: Mean-field approximation 

The interaction term • s° in Eq. Q can be decom¬ 
posed into several quartic operators, which are decou¬ 
pled through a Hartree-Fock approximation (see, e.g., 
Ref.fM)). In what follows, all mean values are taken into 
account to obtain the final mean-field Hamiltonian. 

Using the definition of Si and from Eqs. ([ 2 ]) and 
respectively, the axial component of the interaction term 
becomes 




z(c) 




„z(c)\ 


while the planar component can be written as 

+ (sfX'-” - + (.''“V? + (sr)»U - («.•>(»•*“') 

+ {viK} + - (C/XCc) - (CSU - + iy?,}iy°}, 

with the definitions of (a = 0, x, y, z] {3 = c, /) given by Eqs. ([^ and §. 

Equations (A.l) and (|A.2 ) then lead to 


S, . s? « (sD • S, + (S,) • - (S,) . (s?) + ^(V,/) . V,, -f ^(V,,) . V,/ - i(V,^) . (V,,) 

- \{y^s)v^c - ^{y°c)y°f + \{y^f){y^c)- 


(A.l) 


(A.2) 


(A.3) 


Then, substituting Eq. (A.3) 
Eq. (Q, leads to Eq. ([^. 


into the Hamiltonian, 


In addition, in order to fix the electronic density 
and the number of local magnetic moments, the terms 
-m( c\aC^a - Nric) and e/ ( flf^ - Nrif) must 
also be included in the Hamiltonian, Eq. , with nj = 1. 
These terms represent the constraints which are included 
as Lagrange multipliers, whose values of fi and e/ are de¬ 


termined self-consistently. 

For completeness, we recall that the mean values are 
expressed as 

(Si) = m°(cosQ-Ri,sinQ-Ri,0), (A.4) 

(Si) = —TO°(cosQ-Ri,sinQ-Ri,0), (A.5) 

{K°c) = (y^/) = -yo ( a . 6 ) 
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and where N is the number of lattice sites. Then, the Hamil- 

(V,,) = (V|^) = Fo'(cos(Q-R,),sin(Q.R,),0). (A.7) 

We now perform a discrete Fourier transform on the 
conduction electrons operators (and similarly for the / 
electrons), defined as 


^kcT /—- 


Vn 


y^exp(tk-R^)cicr, 

i 


(A.8) 


JfTi^ 

'Hmf = + ^(^fc+Q ~ A*)Ck+Q^Cj^+Q^ H ^ ^ 

k 

Jm2. 


'^ict'^k+Q^ + H.c.) 


(/kt-^k+Qi + H.C.) + 6/ ^ (/kt-^kt + fl+Qlfk+Ql) 


+ (‘^kt'^kt + '^k+Q;/k+Q; + H.C.) + - ^0 X/ ('^kt-^k+Q; + '^k+Q;/kt + H.C.) 


+ JNrnjm^ + ^JNVq — ^JNVq^ + Nuc^m — NnfCf, 


(A.9) 


where = —2t[cos(fc2,) + cos(/i:y)], while et+Q = —2t[cos(fc2, + q^) + cos{ky + Qy)]. Then, limiting ourselves to 
non-degenerate subspace (k t,k -|- Q ),), where the base vectors are (in a Nambu spinor representation) 




'^k-i-q;’ 


/kf •^k-l-Q^) 


the Hamiltonian is 


'Hmf = E ^|^_^(kt,k+Q|)^^k + const., 

k 


where 


i7(kt,k+Q;) 


/ efe - M 


iJmO 

^k+Q - 

pVo 

\JVi 

\ yvo 

yvo 


yvo 

yy \ 

yvi 

yvo 

o 

-ym° 

-ijm° 

e/ ) 


(A.IO) 


(A.ll) 


(A.12) 


In some instances, the 4x4 Hamiltonian matrix (A.12) 
reduces to simpler 2x2 matrices, saving a significant 
amount of CPU time. One particularly interesting ex¬ 
ample is when the self-consistency process converges to a 
ferromagnetic state, that is, one with Q = (0, 0). In order 
to determine the most general 2x2 Hamiltonian matrix. 


we make use of the rotational symmetry, and take 

(S,) =m5 (0,0,1), (A.13) 

(s^) = -m° (0,0,1), (A.14) 

and 

(W) =Uo'(0,0,1), (A.15) 

on the Hamiltonian of Eg. (|4| . 

Then, taking Eqs. ( ATSfto (A.15) into Eg. (|^, the 
Hamiltonian becomes 
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7m^ 

UmF = -t ^ {c\aCja + H.C.) + ^ ( 


^ E + 7 (3^0 - K) E (4t/.T + Mt) 


+ j(3K) + K) E i4ih + /E 4 ) + - IjNV^ + JNm^fml 


(A.16) 


Fourier transforming the operators in the previous equation, and adding the constraint terms leads to 


aJmPf 


Hmf = E (efc - ^ + Ep^)cLcka + E “ ^^)fljka + 7 E (3^0 - 'tFoO (cL/ka + fL^Ua) 


J ■ 


k.cr 


k,(T 


k,(T 


+ ^JNVq — ^JNVq^ + JNm^jfn^ + fiUcN — ejnjN. 


Using a Nambu spinor representation /j^^), it becomes 

"Hmf = E + const., 


where 


Hika) = 


k(T 




aJmPf 


l{3Vo-aV^ 


;(3Uo-aUoO ef- 


which provides the eigenvalues 


U± =1 

ka 2 

which Cfc = Cfe — /X. 


h + ef + ^(mj-m°) 




aJ 


-[ 2 


h-ef + + m0) 


+ ^(3Uo-aU')' 


(A.17) 


(A.18) 


(A.19) 


(A.20) 


It is interesting to note that the spectra of the 4x4 Hamiltonian are equivalent when Q = (0,0), irrespective 
and of the two 2x2 representations of the mean-field of Vq and Vq' vanishing or not. 
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